Dynamical properties across a quantum phase transition in the Lipkin-Meshkov-Glick 

model 
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It is of high interest, in the context of Adiabatic Quantum Computation, to better understand 
the complex dynamics of a quantum system subject to a time-dependent Hamiltonian, when driven 
across a quantum phase transition. We present here such a study in the Lipkin-Meshkov-Glick 
(LMG) model with one variable parameter. We first display numerical results on the dynamical 
evolution across the LMG quantum phase transition, which clearly shows a pronounced effect of the 
spectral avoided level crossings. We then derive a phenomenological (classical) transition model, 
which already shows some closeness to the numerical results. Finally, we show how a simplified 
quantum transition model can be built which strongly improve the classical approach, and shed 
light on the physical processes involved in the whole LMG quantum evolution. From our results, we 
argue that the commonly used description in term of Landau-Zener transitions is not appropriate 
for our model. 

PACS numbers: 03.67.Lx 



I. INTRODUCTION 

Under a continuous change of the Hamiltonian param- 
eters, a quantum system, initially in its ground state, 
can undergo transitions to excited states. This point 
was already studied in the early days of quantum me- 
chanics, with the celebrated analysis of the two level 
case by Landau[l] and Zener [2]. Here we first display 
some numerical results on the dynamical evolution across 
the Lipkin-Meshkov-Glick (LMG) quantum phase transi- 
tion with two very different pattern whenever the critical 
point is reached from both sides of the quantum phase 
transition (QPT). We then write down a phenomenologi- 
cal (classical) transition model (paragraph HI), which al- 
ready shows some closeness to the numerical results. Fi- 
nally, we show how a simplified quantum transition model 
(paragraph IV) can be built which strongly improve the 
classical approach, and shed light on the physical pro- 
cesses involved in the whole LMG quantum evolution. 
This question is again under a strong focus following the 
raise of the so-called "Adiabatic Quantum Computation" 
(AQC) [3] approach to quantum computation. The basic 
idea here is to initially prepare a quantum system in an 
(easy to prepare) ground state of a known Hamiltonian 
Ho and to encode the answer of a given computational 
problem in the ground state of a final (or "Problem") 
Hamiltonian Hp. The simplest way to implement AQC 
is to take, as time-dependent Hamiltonian H{t)^ a linear 
interpolation between Hq and Hp^ 



H{t) = {l-^)Ho + ^Hp 



(1) 



with t the physical time, and T the total evolution time. 
In the following, we shall use s = t/T ^ with s G [0, 1]. 

According to the quantum adiabatic theorem, if pass- 
ing from the initial to the final Hamiltonian is done slowly 
enough, the system remains in the ground state, and 
therefore, at the end of the evolution, will be in a state 



that encodes the solution of the computational problem. 
A clear distinction between standard (gate-like) quan- 
tum computation (QC) and AQC approaches is the way 
algorithms computational efficiency is quantified. In the 
standard QC case, the complexity refers to the number 
of logical gates to be applied as a function of n^, the 
number of qubits. In AQC, the computational efficiency 
is related to the value that T must be given, such that 
the adiabatic approximation holds, roughly estimated as 



r»A: 



(2) 



where /S^min is the minimum value of the energy gap be- 
tween the ground and the first excited states taken along 
the evolution. The way /S^min (and therefore T) scales as 
function of Uq leads then to the computational efficiency 
in the AQC case. 

Despite these differences, AQC has been proved to be 
(computationally) equivalent to QC [5|, Q. Its main in- 
terest is nevertheless that it provides a direct physical 
implementation, and, to some respect, a different kind of 
robustness to errors and decoherence [7]. 

The QC-AQC equivalence calls for a physical effect 
responsible of AQC failure for highly complex computa- 
tional problems. Since adiabatic quantum evolution is 
strongly sensitive to energy differences from ground to 
excited states, quantum complexity is clearly expected 
when gaps vanish for some value of 5, which is the usual 
signature of a quantum phase transition [8]. Such situ- 
ations with gaps closing exponentially with Uq have in- 
deed been numerically observed for computationally hard 
problems [o, ^. 

One therefore faces the interesting, although expected, 
picture that it is not the ground state itself (for 5 = 1) 
that characterizes the problem complexity, but the na- 
ture of the process, in parameter space (for s < 1). This 
should be related to the already noticed strong relation- 
ship between AQC and Quantum Annealing problems. 



well studied in the past years in the field of complex sys- 
tem (e.g. spin glasses) [10]. 

It is therefore of high interest to study, in the transi- 
tion region, the dynamical and spectral properties in the 
lower part of the spectrum. In particular, one should not 
only focus on the first excited eigenstate coming close 
to the ground state, but in fact to a whole set of ex- 
cited levels. This is in particular expected whenever the 
ground state nature is drastically changed, since its new 
decomposition is mainly weighted by the states for which 
avoided crossings appear along the adiabatic process. So 
the state dynamical evolution has to be thought as a com- 
plex transition cascade, rather than independent Landau- 
Zener (LZ) processes. Therefore, despite its great success 
in other systems, the LZ theory, in particular its ability 
to calculate the transitions between quantum states, may 
not be used in the present AQC case, or be severely cor- 
rected [ll|, [l2[. It should be stressed in addition that LZ 
theory is intrinsically non- adiabatic, which suggest, for 
the AQC slow evolution processes, to go back to more 
standard adiabatic analysis [4]. 

To get a better understanding of these processes, we 
propose here to study the dynamical properties across the 
quantum phase transition in the Lipkin-Meshkov- Click 
model, a simple (solvable) model which exhibits some of 
the expected features of AQC Hamiltonians. The LMC 
model was introduced long ago in 1965 to study shape 
transitions in nuclei [13], and has been, since then, pro- 
posed to describe many systems ranging from interact- 
ing spin systems [14] to Bose-Einstein condensates [15j 
or magnetic molecules such as Mni2 acetate |16j . 

The LMG model describes a set of N spins ^ mutually 
interacting through a XY-\ike Hamiltonian and coupled 
to an external transverse magnetic field h. This Hamil- 
tonian H can thus be expressed in terms of the total spin 
operators Sa = ^i=i ^^/^ where the cr^'s are the Pauli 
matrices: 



l.OP 



H = -j^{l.Sl 



■lySl)-hS, 



(3) 



In the following, we only consider the maximum spin 
sector S = A^/2, with A^ even and A^+ 1 levels. Although 
many different methods have been used to study its ex- 
citation properties, the richness of the full spectrum has 
only be revealed quite recently by means of numerical di- 
agonalizations [17, 18], and then, at the thermodynamic 
limit, in an analytical form [19|. Of interest here is the 
determination of the so-called "exceptional points" in the 
density of states, where the density of states is singular 
and the level separation vanish with N . These points 
gather, as N tends to oo, on a curve which we call here 
the "critical gap curve" (CCC). 

Here, for sake of simplicity, we set 7^ = 0; it is clear 
that, up to a global reparametrization, the Hamiltonian 
only depends on the ratio between 7^^ and h. To express 
this Hamiltonian in the form given by equation ([1]), we 
write 7aj = 5, /i=l — 5, which leads to Hq = —Sz and 
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FIG. 1: Spectrum of the "s" -dependant Lipkin-Meshkov- 
Glick model for N = 20, with 21 levels. 



The energy levels as a function of s are displayed in 
figure [1] for N = 20. At the thermodynamic limit, this 
system undergoes a second order quantum phase tran- 
sition for s = 1/2, whose effect is already visible with 
A/" = 20, in terms of levels pinching. The locus of avoided 
crossing levels appears very close to a straight line (exact 
CCC for infinite N) starting at the QPT for s = 0.5, and 
E/N = —0.5, and reaching £^ = for s = 1. 

In later plots, we shall use the (normalized to one) 
integrated density of states x, in the range [0, 1], instead 
of the energy. The CCC still has a simple expression, 
Xc(s), which reads 



Xc{s) 



1 cot 

TT 






y(i-^)(2^-i) 

(4) 



H, 



- _Q2 



S'JN 



for s e [0.5, 1] 

In the following, we first display some numerical re- 
sults on the dynamical evolution across the LMC quan- 
tum phase transition (paragraph II), with two very dif- 
ferent pattern whenever the critical point is reached from 
one side and another of the QPT. We then write down 
a phenomenological (classical) transition model (para- 
graph HI), which already shows some closeness to the 
numerical results. Finally, we show how a simplified 
quantum transition model (paragraph IV) can be built 
which strongly improves the classical approach, and shed 
light on the physical processes involved in the whole LMG 
quantum evolution 



II. DYNAMICAL EVOLUTION : NUMERICAL 
RESULTS 

A. Forward evolution 

We numerically solve the dynamical (arbitrarily called 
"forward") evolution, proceeding as follows. The initial 
state is the ground state corresponding to the 5 = 



Hamiltonian. The total evolution time T is a multiple 
of a fixed time interval AT, during which the Hamilto- 
nian parameters are kept fixed and the quantum evolu- 
tion is computed by mean of a standard second order 
discretization method. The final state, after AT, serves 
as the initial state for the next step, with (slightly) var- 
ied Hamiltonian parameters. Therefore, the larger the 
T value, the smaller the effective Hamiltonian variations 
from one step to the next, and therefore the closer to an 
adiabatic evolution. Another parameter is A/", the sys- 
tem size. Increasing TV decreases the gaps, which eases 
the transitions to excited states. Computations are done 
here with N = 50 and three T values, corresponding 
qualitatively to fast (T = 1), medium (T = 50) and slow 
(T = 100) evolutions. 

The levels occupancy, as a function of s, are displayed 
in Figure [2] (upper plots). Also shown is the CGC curve, 
in order to track the role of the gap closing phenomenon 
in the quantum evolution. 

As can be clearly seen, a common feature of these evo- 
lutions is that the system almost remains in its ground 
state before reaching the quantum phase transition re- 
gion. 

Then, not only do the quicker evolutions drive the sys- 
tem to excited states transitions, but this evolution is 
clearly controlled by the position of the avoided cross- 
ings, as marked by the critical gap curve. 

As expected, for slower evolutions (larger T), the 
ground state is not completely depleted, its population 
oscillates with time (as seen on the figure) and eventually 
stabilizes (see for example Ref. [llj). 



which is first visited in the backward evolution, the gaps 
are smaller with respect to those in the left part. As a 
consequence, for the same T value, transitions to excited 
states are more probable, and start occurring even before 
crossing the quantum phase transition (Figure [2] (lower 
left)). 

Note finally (Figure [2|) that the states population is 
slightly less peaked (it displays a sudden small anti- 
bump) when the CGC curve is crossed. This is again 
a manifestation of the enhanced transition probability 
along this curve. 



III. A PHENOMENOLOGICAL MODEL FOR 
THE QUANTUM EVOLUTION 

We shall be interested, in this chapter and the next 
one, in computing approximate values for the evolution 
of the probability distribution among the different states. 

We call {|^i(s))} the instantaneous eigenbasis of the s- 
dependant Hamiltonian, H{s)\r]i{s)) = ei{s)\r]i){s))^ and 
write the current state|<l>(5)), in this basis : |*^(s)) = 
'^iCii{s)\rji{s)). We aim to compute the probability 
Pi{s) = \{r]i{s)\^{s))\'^ for the system to be in the i-th 
instantaneous eigenstate at time s. We shall first sup- 
pose that transitions only occur, at the same rate, from 
an instantaneous eigenstate \r]i) toward states |/?i+i) and 
\{r]i-i). The transition rate matrix F inherits a tridiag- 
onal form, but with nevertheless 5— dependant elements. 
This leads to the following differential equation for the 
probability 



B. Backward Evolution 

The above observation that the whole spectrum infiu- 
ences the overall quantum state evolution, together with 
the fact that the LMG spectrum is far from being sym- 
metrical (see figure [1]), leads to expect a qualitatively 
different time evolution whenever the system is driven 
backward, which reads 



Hinv{t) = sHo^{l-s)Hp 



(5) 



The levels occupancy, displayed in Figure [2] (lower 
plots), indeed shows a very different pattern. A first ex- 
planation arises quite naturally : in the forward case, the 
system encounters the minimal energy gaps in an ordered 
sequence that allows the current wave function to spread 
in the spectrum, with a high probability to change its 
eigenstate decomposition along the avoided crossings. In 
the latter (backward) case, once the system encounters 
the first small gap, and possibly leaves the ground state, 
it never meets again avoided crossings situations, and 
therefore do not proceed significantly to higher energies. 

In addition, the levels population displays other quali- 
tative features which can be explained by looking to the 
spectrum. In the right part of the spectrum (Figure [1]), 



P^{s) = r,+i,,(5)p,+i(5) + r,_i,,(5)p,_i(5) 

-(r,,,+i(5) + r,,,_i(5))p,(5) (6) 

We choose a form for the rates Ti^i±i{s) which follows 
a generic adiabatic prescription 



r«,j(5) = 



Tb 



(A^'^isW 



(7) 



where A*'-^(s) is the instantaneous gap between levels i 
and j, calculated from the spectrum, T is the evolution 
time and b a (adjustable) coupling parameter . 

Numerical solutions of the master equation (|6]) are dis- 
played in Fig. [3] (up left). This (very) simple model fairly 
reproduces some features of the computed evolutions. In- 
deed, for the forward evolution, we find the sequential 
transition driven by the CGC curve while for the back- 
ward evolution, a sudden transition to excited states and 
a subsequent saturation effect (after crossing the CGC) 
are recovered. 

However, in the backward case for instance [ Fig. [H - 
(up left) ], this simple model fails in describing correctly 
the lowest levels occupations. Indeed, the ground state 
remains here the most populated state during the evo- 
lution, a feature which is clearly not found in the full 
numerical simulations. 
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FIG. 2: Forward and backward evolution of the level populations for a system of size A/" = 50 for different the total evolution 
times T. The red line is the analytical CGC curve(see text). The three upper figures correspond to forward evolution : left, fast 
evolution (T = 1); middle, intermediate speed (T = 50); right, slower evolution (T = 100). The three lower figures correspond 
to the backward evolution, with the same speed (T values) as the plots above 



At this point an important remark must be done. 
In this phenomenological model, we choose a transition 
rate ([7j) which depends on the inverse of the square 
of the energy gaps. This choice of an "adiabatic"-like 
transition rate is, in a certain sense, arbitrary. An- 
other possible choice would have been to take Tij{s) ex 
exp {— 6(A*'-^(5))^} (with b > 0), which mimics the 
Landau-Zener transition rate. We tried to plug this type 
of behavior in the rate equation, but could not find any 
reasonable agreement with the numerical results. This 
suggests that, for the present model, a generic Landau- 
Zener-like description is not appropriate. 



IV. SIMPLIFIED QUANTUM MODEL FOR THE 
ADIABATIC EVOLUTION 



In order to better describe the full dynamical process, 
we need to improve the previous approach, and incorpo- 
rate quantum effects more precisely, as follows. If the 
evolving quantum state is written in the instantaneous 
basis of H{s)^ with the gauge choice {Tjn{s)\T]n{s)) = 0, 
the am{s) coefficients satisfy the equation 



^ f ^ ■ f \ f \ \^ {r]m{s)\dsH{s)\r]n{s)) . . 

dsam[S) = -iem[S)am[s)- 2_^ An,.rr^.r.\ ^n[S) 



n^m 



A"'"'(s) 



(8) 

where, as above, A"^'^(s) is the instantaneous gap be- 
tween levels n and m. Note that solving equation ([8]) 
requires the knowledge of both the time-dependant eigen- 
values and the matrix elements {irim{s)\dsH{s)\irin{s)). 

Let us further simplify this equation by using the nu- 
merically derived gaps and use an approximate form for 
the matrix elements. In addition, we limit the transi- 
tions from level n to neighboring levels n + 1 and level 
n — 1, leading to an effective ID quantum chain model. 
As can be seen in figure [5l this is justified far from the 
CGC curve because matrix elements clearly display an 
exponentially decay form with distance in the effective 
chain. 

We then use two forms for the matrix elements. In 
the simplest case, we take them equal for any inter- level 
transition and constant with s. This already improves 
the previous phenomenological approach, as can be seen 
in figure [3] (up-right). To further improve our model, 
we numerically compute the matrix elements in equa- 
tion ([8]). We find that in general, for a given level n, 
{'nn{s)\dsH{s)\rin+i{s)) have a maximum near the CGC 
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FIG. 3: Levels population in the forward evolution case for a 20 level system and T = 10. Phenomenological rate equation 
{T h — 0.01) (up left), effective quantum chain with constant matrix element (up right), effective chain with improved matrix 
elements (down left), full numerical evolution (down right). 



curve and then show a fast decrease. We fit this 5- 
dependence with a Gaussian form, centered at the crit- 
ical point. In addition we find that these matrix ele- 
ments have maximal values, as a function of n, which 
can be well approximated by a logarithmic behavior. 
We therefore write them as {irin{s)\dsH{s)\irinJri{s)) = 
(a -|- 61ogn) exp (—7(5 — 5o(n))^), where a, b and 7 are 
fitted parameters and so(n), the s value where the gap 
between levels n and n + 1 is minimum, is very closely 
approximated from the analytic expression of the CGC 
curve. 

The numerical simulation with this latter approxima- 
tion is shown in figure [3] (down left). It shows a clear 
improvement with respect to the phenomenological ap- 
proach ([6]), and even to the above constant matrix el- 
ement approximation, in particular after crossing the 
CGC line, where the depleted population of the lowest 
levels is better reproduced. 

Backward evolution, treated with the same approxi- 
mations, are presented in figure HJ with similar trends as 
in the forward case. 

Let us stress that the qualitative form of the computed 
evolution does not depends critically on the s dependence 
of the matrix elements. Analogous results are obtained 
within a vast range of 7 coefficients and even with a differ- 
ent functional dependence, as long as the approximating 
function remains well peaked around the critical point. 



This confirms that the CGC drives the main feature of 
the dynamics; on the other hand, the absolute value of 
the matrix elements determine the "fine details" of the 
evolution, such as the ratio of the population levels. 

A final remark concerns one important basic assump- 
tion of the above approaches (both phenomenological and 
quantum), that this system is well approximated by an 
effective chain with only nearest neighbor transitions be- 
tween levels. This is true only far from the CGC curve; 
near the critical point, long "distance" transitions occur. 
This point has been analyzed in a semi-classical frame- 
work, and will be presented elsewhere [20]. 



V. CONCLUSION 

We have studied time dependant dynamics of the 
Lipkin-Meshkov- Click model driven across its Quantum 
Phase Transition Point. The dynamics of the quantum 
evolution, not restricted to the lowest level occupancy, is 
determined by the spectral critical gap curve, where the 
energy gaps vanish at the thermodynamic limit. In order 
to compare with the full numerical solution, we have de- 
veloped simplified models for the transitions during the 
evolution. 

First, we use a phenomenological rate equation ap- 
proach, with adiabatic-like transition rates, which al- 
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FIG. 4: Levels population in the backward evolution case for a 20 level system and T = 10. Phenomenological rate equation 
{T h — 0.01) (upper left), effective quantum chain with constant matrix element (upper right), effective chain with improved 
matrix elements (lower left), full numerical evolution (lower right). 




0.001 



FIG. 5: Logarithmic plot oi M(k) = \{rin+k{s)\dsH{s)\rin{s))\ 
(with n = 10 and A/" = 20 ) as a function of the levels distance 
k and for different s (from lower to upper curves: s = 0.75, 
s = 0.85 and s = 0.90), together with linear fits. 



ready recovers the role of the CGC curve in driving the 
main quantum evolution. But this approach misses some 
important features of the quantum evolution. 

We then improve our description by building a quan- 
tum model, which treats the inter- level jumps as the con- 
sequence of an effective interaction between the instanta- 
neous levels of the "s" -dependant Hamiltonian, restricted 



to nearest-neighbor level interactions. In a first step, this 
interaction is only varied following the values of the "s"- 
dependant gaps, which already compares better with the 
full numerical solution. We then further improve this 
effective chain model by including an approximate form 
for the rate of change of the Hamiltonian averaged over 
neighboring levels. In that case, the main features of the 
quantum evolution are recovered. 

Future investigations should focus on a finer descrip- 
tion of the quantum evolution near the critical curve. In 
particular, long range interactions between instantaneous 
levels, which come close in energy near that curve, have 
to be taken into account. 
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Note added - While finishing this manuscript we 
became aware of a recent work by Caneva et al. ^l| 
where the adiabatic dynamics for the LMG model is 
also studied. They focus on computing and analyzing 



the residual annealing energy as a function of the 
total adiabatic time. In particular, they numerically 
identify three different regimes for the system evolution. 
For intermediate adiabatic times, the Kibble- Zurek 
mechanism, based on an analysis of "defects" and 
domain walls, used in [22] for one-dimensional systems is 



invoked to analyze the computed behavior. We may just 
remark here that this may be questionable : the LMG 
model corresponds to a fully connected set of spins 1/2, 
therefore asymptotically in infinite dimension, and it is 
not clear how to adapt the Kibble- Zurek approach in 
that case. 
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